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The new version Q of the HERWIG event generator Q allows one to simulate 
hadronic final states in electron-positron annihilations using the O(a^) partonic 
matrix elements (MEs) for e~^e~ — > qqgg and e~^e~ — > qqQQ. The consequent phe- 



> . 

, nomenology is discussed in the case of four-jet signatures of particular relevance to 

future electron-positron Linear Colliders (LCs), such as in Higgs boson searches. 
OO I The new approach adopted in HERWIG is compared to the standard 0{as) de- 

' scription, as well as to other recently presented event generators also based on 

■ e~^e~ 4 parton MEs. 



Ph| Four-jet events will play a crucial role at future LCs 0. Just one example is sufficient to 

Q^l make this point. Such machines are an ideal laboratory where to pursue studies in the 

Higgs sector [Q. In this respect, it should be recalled that e~^e~ 4 jet processes from 
QCD will represent a serious noise in the search for Higgs particles, both in the standard 
electroweak (EW) theory and in possible extensions. In particular, in the Standard Model 



^ I (SM), the dominant Higgs production channel proceeds via the 'bremsstrahlung' process 

e'^e~ HZ P], followed by the hadronic decay H —>■ bb (in the intermediate regime), 
with the Z boson going into pairs of jets about 70% of the times. 

Thus, it is very important that four-jet events from QCD are correctly implemented 
in the Monte Carlo (MC) programs that are widely used in phenomenological studies 
of hadron production at e~^e~ colliders. In this connection, it has always been rather 
worrying that certain aspects of four-jet production were never well described by the 
standard ^0{as) ME -|- parton shower (PS)' MC programs. As an illustration, one may 
recall the four-jet studies performed by the ALEPH Collaboration at LEPl 0, which had 
revealed a significant disagreement between data and MCs in the typical four-jet angular 
variables (see Ref. ^ for the definition and the discussion of some of their properties): (i) 
the Bengtsson-Zerwas angle xbz; (ii) the Korner-Schierholz-Willrodt angle $ksw; (iii) the 
(modified) Nachtmann-Reiter angle ^^rj i^^) the angle between the two least energetic 
jets 6134. 

^Talk given at the 2nd ECFA/DESY Study on Physics and Detectors for a Linear Electron-Positron 
Collider, Padova, Italy, 4-8 May 2000. 
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Figure 1: Differential distributions for the polar angle of the reconstructed Z. 



In Ref. it was argued that the discrepancy was nothing but the evidence that 
the standard ^0{as) ME + PS' MCs do not provide a correct description of the spin 
correlations among the various partons, particularly at large jet separations. In fact, the 
above angular quantities are defined in terms of the three-momenta of the particles in 
the final state and are precisely the spins of the latter that induce very peculiar orien- 
tations in (i)-(iv). Conversely, ME' programs (for example, the 'O(a^) ME + 
string fragmentation model' as implemented in JETSET, i.e. with no parton shower) 
yielded a much better angular description of four-jet final states, see Refs. In fact, 
all spin correlations are naturally included in a full Ola^) ME calculation, but they are 
not necessarily present in a PS emulation based on the 0{as) scattering, although their 
availability (in the infrared limit, where soft and coUinear correlations can be factorised in 
analytic form) is a feature of some of the above mentioned MC codes (for example, for the 
HERWIG implementation, see Ref. |iy]). However, even an '(9(q;^) ME + fragmentation' 
model, without the PS evolution, is inadequate to describe QCD four-jet production in 
e'^e~ collisions. The problem is that such ME models contain 'ad-hoc' hadronisation 
which is adjusted to produce a good agreement with some data [^, |TT|, but they cannot 
be extrapolated to other energies. Furthermore, their description of the sub-jet structure 
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Figure 2: Differential distributions for the opening angle of the jets in the di-jet system recon- 
structing the Z. 



is very poor (see Ref. [12] for an overview 



That such deficiencies in the description of four-jet final states could be cured by an 
'^(a^) ME + PS' (plus hadronisation) approach has been clear for some time. Now, 
several event generators have been made available, that resort to this implementation: 
HERWIG 6.1 111], PYTHIA 6.1 [|l|], 4JPHACT and APACIC++ It is the intent 
of this note to describe the impact that the new description of four-jet final states can 
have in some contexts which are particularly crucial to the physics of future LCs. To do 
so, we will borrow as an example the case of the HERWIG event generator, though many 
of its features can equally be appreciated in the case of the other MCs. Eventually, we 
will also highlight the differences and similarities among these various programs. 



Ref. |]T6[ has already made the point, for the case of LEP, that systematic differences 
exist between the use of the 0{as) and 0{al) MEs in the simulation of QCD four-jet 
events. This is true not only for the case of QCD studies (e.g., in the determination of the 
QCD colour factors) but also in the context of Mw± measurements. Here, we illustrate 
similar effects taking place in a higher energy environment, such as a LC operating at 

= 500 GeV. 
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Figure 3: Differential distributions for the invariant mass recoiling against the reconstructed Z. 



Before proceeding to the numerical part of our note, we would like to spend a few 
words about some technical aspects of the HERWIG implementation. The new pro- 
cess describing electron-positron annihilation into four jets can be simulated by setting 
IPR0C=600+IQ, where a non-zero value for IQ guarantees production of quark flavour IQ 
whereas IQ=0 corresponds to the natural flavour mix. (IPR0C=650+IQ is as above but 
without those terms in the MEs which orient the event with respect to (w.r.t) the lepton 
beam direction.) The MEs are based on those of Ref. []T^| with orientation terms from 



ISf . Infrared (soft and collinear) divergences are avoided by imposing a minimum sepa- 



ration among the four partons (hereafter, denoted by ?/inf - Y4JT in the code, by default 
equal to 0.01). At present, the jet separation is calculated using either the Durham or 
Jade metrics (see, e.g. Ref. for deflnitions and a review on jet clustering algorithms). 
This choice is governed by the logical variable DURHAM (whose default is .TRUE.). Note 
that the phase space is for massless partons, as are the MEs, though a mass threshold 
cut is applied. The scale EMSCA for the parton showers is set equal to ^/sy^, where y^iri 
is the least distance, according to the selected metric, between any two partons. Finally, 
the jet resolution y^ut ^ used to select a multi-jet flnal state - should be chosen such that 
it is always larger than ?/inf. 



3 



An ambiguity arises in the treatment of interference terms between diagrams with 
different colour structures. In fact, apart from the trivial case e+e^ QqQQ, with 
q Q, the MEs for e~^e^ QQgg and e~^e^ — > qqqq (i.e. q = Q) receive contributions 
from two distinct colour flows each. In general, the interference terms between the latter 
are not positive definite and a choice has to be made, on how to assign them to one or 
the other of the two leading colour contributions (although suppressed by w.r.t. 
the latter, they cannot be neglected without inducing sizable effects in both total and 
differential rates). For illustration, we describe here the treatment adopted in HERWIG 
for e^e~ — > qqgg and refer the reader to the forthcoming new HERWIG manual (see 
also [21]) for details on the e^e~ — > qqqq subprocess. Things go as follows. Two-quark- 
two-gluon final states are produced via the following eight Feynman graphs (here, 5 labels 
the e"'"e~ — >• 7, Z current): 



Ti 




To 





T4 








The two aforementioned colour flows are identified by the terms proportional to 
{t'^t^)i^i^ and {t^t^)i^i^, where A and B are the colours of the gluons labelled 1 and 
2, respectively. In fact, it should be recalled that the colour pieces associated with the 
triple-gluon vertices are nothing but the structure constants /^-^^ of the SU{Nc) gauge 
group, for which 

[t it \i3i4, = {t t — {t t = if ^4314! (2) 

X being the label of the intermediate gluon and isyii the outgoing quark colour indices. 
Therefore, one can conveniently group the original eight Feynman amplitudes as 

^+ = E^^ + E^- M_ = ^T,-^T,. (3) 

i=l 1=7 1=4 1=7 
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The two positive definite contributions and corresponding to the two fundamental 
colour connections needed for the interface to the HERWIG parton shower (see Ref. 
can, for example, be obtained as: 

N+= D (|M+|2-i^|M+ + M_ 

c 



|2 


|M+ 


2 




P + 


|M_|2 


|2 


|M_ 


2 




l' + 


|M_|2 



where the colour pre-factor is (A^c" = 3) 

D = ^{Nl-ll (5) 
so that the total squared amplitude can be recovered from 

M\qqgg) = Y,N,. (6) 
i=± 

In the right-hand sides of eq. (^, the first terms are the 'planar' whereas the second 
are the 'non-planar' ones Here, would correspond to the t-channel colour fiow 



(4 2, 2 1, 1 ^ 3) and A^. to the w-channel one (4 1, 1 ^ 2, 2 ^ 3), that is, 2413 



and 4123 in the notation of Ref. This is the default procedure adopted in HERWIG 
6.1 that we have maintained in our analysis (see Ref. [||] for other possible choices). 

Several kinematic quantities can be used to disentangle Higgs events from the back- 
grounds at a LC. One example is the polar angle of the reconstructed Z boson. In 
fact, in the case of the signal, one expects the gauge vector to be very central, own- 
ing to the s-channel and spin dynamics of the Higgs process 0. In contrast, no such 
a behaviour would appear in the QCD noise. Fig. |I] shows this variable as obtained by 
using 0{as) (from HERWIG 5.9, IPR0C=105) and 0{al) (from HERWIG 6.1, IPR0C=605) 
MEs, supplemented by the subsequent PS and the hadronisation stage. The four-jet final 
state has been selected by using the Durham (D) jet-clustering algorithm with resolution 
Vent = 0.0010. Furthermore, we have imposed that one pair of jets reproduce the Z mass 
within 10 GeV, i.e. |Mjj — Mz\ < 10 GeV, as illustration of a typical selection proce- 
dure of e'^e" HZ — > 4 jet events. These are the two jets that we will consider to be 
produced in the Z ^ 2 jet decay. For reference, we also have superimposed the same 
angular spectrum as obtained from the Higgs process, again generated by using HERWIG 
(IPR0C=305). As self-evident from the choice of the IPROC numbers, we always require 
two 6-jets in the final state, for both QCD and Higgs processes. The HERWIG default 
choice Mh = 150 GeV is used for the SM Higgs mass. Also notice that we normalise the 



^Notice that exactly four jets are required to be reconstructed at hadron level. To allow for n > 4 jets 
in the final state (eventually clustered into four) would not alter our conclusions. Note that we use the 
so-called E recombination scheme: two tracks are clustered by summing their four-momenta. 
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y(cut) Durham-E at Parton Level 



• Herwig 6.1, starting by 4-parton ME 
o Pytliia 6.1, starting by 4-parton ME 
A 4jphact + Pythia 6.1 
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y(cut) Durham-E after Hadronization 



FigTirc 4: The Ucut dependence of foTir-jct events at parton level before (top) and after (middle) 
showering effects as well as at hadron level (bottom). 



cross sections of all processes to unity. In other terms, we implicitly assume that the two 
QCD descriptions produce the same event rate, for a given choice of qjs, jet- clustering 
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scheme and ?/cut- In reality, there exist differences in the total production rates between 
the two QCD implementations, that could be source of further systematic uncertainties. 
However, we leave these aside for the time being, as we only concentrate on the kinematic 
effects. It is clear from Fig. |l| that not only the two QCD descriptions show significant 
differences between each other in the polar angle distribution of the di-jet system emu- 
lating a Z decay, but these also persist where the Higgs process accumulates. Similar 
effects can be appreciated in the distribution of the relative angle between the two jets 
reconstructing the gauge vector, see Fig. |. This quantity too is used in the selection of 
Higgs events at a LC ||^ . Thus, although the behaviour of the two QCD implementations 
is rather similar in the vicinity of the Higgs peak, see Fig. ^, the effect of cuts enforced on 
cos^bcam aud/or cos^jj can be dramatically different for the actual number of background 
events falling around Mj-ecoii ~ ^u- This can be appreciated in Tab. [I|. 



HZ 






72. % 


9.8 % 


12. % 



Table 1: Percentage of four-jet events in a window |Mrecoii — Mh\ < 30 GeV, after the 
cuts I cos6'bcam| < 0.6 and cos6'jj < 0.5. 

As mentioned earlier on, at least three other MCs exist that have the option of gen- 
erating multi-jet final states starting from the 0{a'^) MEs. Among these, we consider 
here PYTHIA 6.1 and 4JPHACT. The former makes use of the same (massless) MEs 
implemented in HERWIG (with a similar approximate procedure to account for mass 
effects fl^), the main differences being in the treatment of the PS and hadronisation. 



The latter, while using showering and fragmentation dynamics borrowed from PYTHIA, 
employs the fully massive MEs of Ref. p3[ and a massive four-body phase space as well. 
A comparison among these MCs has been performed by Silvia Bravo for ALEPH and we 
present here just a selection of the results she obtained^. This is done in Fig. ^, where the 
I/cut dependence of the four-jet rates (all flavour combinations), as obtained in the three 
MC environments, is presented at three different stages. First, at the four-parton level 
before showering effects, then after the latter have taken place, finally after hadronisa- 
tion. The comparison is performed at y/s = Mz and exactly four jets are required to be 
selected by the Durham jet finder (using E-scheme recombination), for a given ?/cut value. 
Normalisation is to the number of generated events: 1,000,000 in each case. The overall 
agreement is quite good, particularly at the experimentally more relevant hadron level. 
In the end, systematic errors induced by the different treatment of MEs, phase space, PS 
and hadronisation in the three MCs are typically much smaller than those between any 
0{a'^) description and the standard 0{as) ones. 



Others can be found at http://www.pd.infn.it/ecfa/. 
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In summary, we have shown that there exist sizable differences in the MC simulation 
of key LC phenomenology in four-jet final states, depending upon whether e~^e~ qqg 
or e+e^ 1<i9g plus e"'"e~ — » qqQQ partonic scatterings are used to initiate the PS of 
the event generator. However, several reliable MC programs exist by now, each exploiting 
a different implementation of either approach. Thus, new, more sophisticated analyses 
of four-jet events will soon be available, these allowing one to pin down to even better 
accuracy the physic potential of a future e+e~ collider operating at the near-TeV scale. 

Before concluding, we would like to briefly remark on other aspects that enter the 
description of hadronic final states with arbitrary jet multiplicity, that were beyond the 
scope of this note. Here, we have been concerned with what we should like to call the 
procedure of 'interfacing' MEs to the subsequent PS, that is, of generating hadronic events 
through higher-order MEs, but limiting the event generation to a phase space region 
where the former are reliable, without any attempt to resort to the latter otherwise. 
In fact, another important issue is that of 'matching' MEs and PS, by which we mean 
the procedure of exploiting a combined approach that uses MEs for large phase space 
separations, PS in the infrared limit and modified MEs (by Sudakov form factors) in the 



overlapping region [Q. This is what is currently done within the APACIC++ event 
generator |15[, which uses exact leading-order (LO) MEs for e^e^ n partons, with 
n < 50. Note that the procedure of Ref. [0] is correct to next-to-leading logarithmic 
(NLL) accuracy and could possibly be extended to next-to-leading order (NLO) as well 
PB| . The interface/matching to the PS of higher-order MEs including loops [^, p7| 
(hence, generating negative weights) remains instead a challenge for the future ||28|| . 
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